The  optical  properties  of  aerosols 


Final  Technical  Report 
by 

F.  Borghese 
September  1997 


United  States  Army 

EUROPEAN  RESEARCH  OFHCE  OF  THE  U.  S.  ARMY 
London,  England 


Contract  number  N68171-96-C-9056 
Contractor:  Centro  Siciliano  per  le  Ricerche  Atmosferiche 
e  di  Fisisca  dell'Ambiente 


Approved  for  public  release;  distribution  unlimited. 


AppsToml  ta  ’  I 


19911112  013 


The  optical  properties  of  aerosols 
Final  Technical  Report  on  Contract  N68171-96-C-9056 
September  1997 


The  research  performed  under  the  present  contract  has  dealt  with  the  follow¬ 
ing  items: 

1)  Optical  properties  of  inclusion-containing  hemispheres  deposited  on  a  per¬ 
fectly  reflecting  surface. 

2)  Optical  resonances  of  homogeneous  spheres  containing  an  eccentric  spherical 
inclusion . 

3)  Optical  properties  of  aggregated  spheres  deposited  on  a  dielectric  surface. 

These  researches  are  fully  expounded  in  papers  that  have  already  been 
submuitted  for  publication  or  are  in  an  advanced  stage  of  completion.  There¬ 
fore,  in  the  present  report,  we  outline  the  motivations  behind  each  work  and 
summarize  the  main  results. 

1)  Optical  properties  of  inclusion-containing  hemispheres  deposited  on  a 
perfectly  reflecting  surface.  The  aim  of  this  research  is  the  assessment  of  the 
extent  to  which  the  presence  of  inclusions  modify  the  scattering  pattern  from 
otherwise  homogeneour  hemispheres  deposited  on  a  perfectly  reflecting  sur¬ 
face.  The  practical  interest  of  this  kind  of  work  is  easily  understood  when  one 
thinks  of  the  hemispheres  as  liquid  droplets  deposited  on  a  metallic  surface  and 
containing  small  pollutant  particles.  In  this  respect  let  us  recall  that  a  perfectly 
reflecting  surface  is  the  widely  accepted  model  for  a  metal  surface. 

In  performing  our  task  we  took  advantage  of  our  previous  experience  in  the 
study  of  the  optical  properties,  and  in  particular  of  the  optical  resonances,  of 
single  and  aggregated  homogeneous  hemispheres  deposited  on  a  metallic  sur¬ 
face.  Indeed  the  technique  that  we  used,  both  in  our  previous  and  in  the  present 
research  is  based  on  Image  Theory,  that,  for  the  case  of  a  perfectly  reflecting 
surface,  yields  particularly  simple  and  effective  computational  procedures.  In 
fact,  we  were  able  to  calculate,  without  undue  computattional  effort,  the  pattern 
of  the  scattered  intensity  as  a  function  of  the  angle  of  observation  for  several 
choices  of  the  angle  of  incidence  and  of  the  polarization  of  the  incident  wave. 
Several  morphologies  were  considered  such  as  hemispherical  and  spherical  in¬ 
clusions  as  well  as  aggregations  of  two  identical  spheres.  The  polarization- 
dependent  patterns  proved  useful  to  get  information  on  the  possible  anysotro- 
py  of  the  inclusions ,  whereas  the  patterns  obtained  for  different  choices  of  the 
angle  of  incidence  give  information  on  the  position  of  the  inclusions  within  the 
host  hemispheres. 

Ultimately,  the  calculated  patterns  give  enough  information  as  to  allow  us  to 
draw  reliable  conclusions  both  on  the  morphology  and  on  the  position  of  the  in¬ 
clusions.  Therefore  we  plan  to  extend  this  work  to  include  averages  both  over 
the  position  and  the  orientation  of  the  inclusions  in  order  to  get  information 
that  may  be  of  more  direct  use  to  the  ineterested  experimentalists. 


2)  Optical  resonances  of  homogeneous  spheres  containing  an  eccentric  spher¬ 
ical  inclusion .  The  study  of  the  optical  resonances  has  proven  to  be  a  reliable 
method  to  get  information  on  the  internal  structure  of  small  particles.  Neverthe¬ 
less,  the  interpretation  of  the  experimental  data  can  be  based  on  general 
grounds  only  for  spheres  either  homogeneous  or  radially  non  homogeneous. 

For  such  objects,  indeed,  as  a  consequence  of  the  spherical  symmetry,  the  reso¬ 
nances  do  not  depend  either  on  the  direction  of  incidence  or  on  the  polarization 
of  the  incident  light.  For  non  spherical  particles  the  resonance  spectra  are  rather 
more  complicated,  as  the  lack  of  spherical  symmetry,  implies  a  dependence  of 
the  resonance  peaks  both  on  the  direction  of  incidence  and  on  the  polarization. 
Up  today  no  general  theory  exists  on  which  the  interpretation  of  the  observed 
spectra  can  be  reliably  based.  For  this  reason  a  theoretical  study  of  the  behavior 
of  the  resonances  of  objects  lacking  the  spherical  symmetry  may  be  of  interest. 
Our  research  on  the  resonances  deals  with  the  calculation  of  the  extinction  spec¬ 
trum  from  homogeneous  spheres  containing  a  spherical  inclusion  that  is  al¬ 
lowed  to  move  from  a  concentric  to  an  eccentric  position.  This  case  is  of  interst 
because  the  calculated  spectra  permit  to  follow  the  behavior  of  the  resonance 
peaks  when  the  symmetry  of  the  scatterer  as  a  whole  lowers  from  spherical  to 
cylindrical.  Our  actual  calculations  were  performed  on  homogeneous  spheres 
of  MgO  containing  an  inclusion  of  solid  Mg:  in  other  words  we  considered  a 
metallic  sphere  (possibly  eccentrically)  coated  by  its  oxyde. 

As  expected  on  the  ground  of  Group  Theory  we  found  that  the  resonances  un¬ 
dergo  a  splitting  when  the  symmetry  lowers.  The  splitting  persists  even  when 
one  considers  the  resonances  from  a  dispersion  of  randomly  oriented  scatterers. 
When  the  latter  spectrum  is  compared  with  those  for  a  dispersion  of  likely 
oriented  objectswe  found  that  it  is  possible  to  get  information  on  the  fraction  of 
scatterers  that  are  likely  oriented.  It  is  clear  that  such  a  possibility  may  be  of 
help  to  study  those  processes  that  imply  the  orientation  of  anisotropic  particles 
by  electric  or  magnetic  fields. 

3)  Optical  properties  of  aggregated  spheres  deposited  on  a  dielectric  surface. 
Image  theory  proved  in  time  to  be  the  method  of  choicefor  the  description  of 
the  optical  properties  of  particles  in  the  vicinity  of  a  perfectly  reflecting  surface. 
Nevertheless,  when  the  surface  is  a  dielectric  one,  i.  e.  when  it  separates  two 
different  media,  the  choice  of  a  different  approach  is  in  order.  Indeed,  several 
papers  showed  that,  in  the  latter  case,  image  theory  yields  a  complicated  proce¬ 
dure  unless  suitable  approximations  or  assumptions  on  the  size  of  the  particles 
are  made.  Our  Group,  working  under  the  preceding  Contract  was  able  to  for¬ 
mulate  a  general  theory  for  the  reflection  of  vector  multipole  fields  of  arbitrary 
multiplicity,  on  a  dielectric  surface.  This  theory  has  then  been  applied  to  the  de¬ 
scription  of  the  scattering  properties  of  spheres,  not  necessary  small,  in  the  vi¬ 
cinity  of  a  dielectric  surface.  The  results  obtained  through  our  approach  proved 
to  be  in  excellent  agreement  with  the  results  of  ab  initio  simulations  and  with 
the  experimental  data  so  that  we  were  encouraged  to  extend  our  method  to  de¬ 
scribe  the  scattering  properties  of  aggregated  spheres.  In  this  sense  our  present 


research  is  a  continuation  of  the  work  performed  under  the  preceding  Con¬ 
tract.  The  motivations  that  led  us  to  perform  this  kind  of  extension  are  easily 
understandable.  It  is  well  known,  indeed,  that  the  solid  particles  in  many  aero¬ 
sols  may  aggregate  to  form  bigger  particles  that  tend  to  fall  under  the  effect  of 
gravity  and  to  deposit  on  any  suitable  surface.  A  reliable  computational  meth¬ 
od  that  allows  to  foresee  the  changes  undergone  by  a  clean  dielectric  surface 
when  aggregated  particles  are  deposited  is  thus  of  interest. 

The  formulation  of  our  theory  is  such  as  to  permit  to  perform  orientational  av¬ 
erages  over  the  particles  and  thus  to  describe  the  situation  that  is  likely  to  occur 
in  actual  experiments.  Nevertheless,  due  to  the  mathematical  complications 
arising  from  the  need  to  describe  aggregated  spheres,  our  calculations  are  not 
yet  complete.  We  can  state,  however,  that  the  preliminary  result  show  the  cor¬ 
rect  limiting  behavior  when  the  surface  tends  to  become  perfectly  reflecting  as 
well  as  when  the  interparticle  interaction  tend  to  weaken  with  increasing  inter¬ 
particle  distance. 

To  complete  this  Report  we  include  the  preprint  of  the  paper  "Optical  resonanc¬ 
es  of  spheres  containing  an  eccentric  spherical  inclusion"  that  has  been  submit¬ 
ted  to  Journal  of  Optics  as  well  as  reprints  of  the  papers  that  have  been  pub¬ 
lished  in  the  present  year  although  performed  under  the  preceding  Contract. 
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Abstract.  The  resonance  spectrum  of  small  spheres  of  Mg  with  a  spherical  coating  of 
MgO  is  calculated  in  the  framework  of  the  transition  matrix  approach.  The  spectrum 
for  concentric  coating  displays  a  single  resonance  peak  whereas  when  the  coating  is 
eccentric  the  spectrum  complicates  as  the  number  of  resonances  increases.  In  the  latter 
case,  the  spectrum  undergoes  evident  qualitative  changes  according  to  the  choice  of 
the  polarization.  Comparison  of  the  spectra  for  a  dispersion  of  particles  oriented  alike 
with  the  one  for  randomly  oriented  particles  suggests  a  possible  procedure  to  evaluate 
the  fraction  of  particles  that  are  likely  oriented. 
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1.  Introduction 

The  optical  resonances  of  spherical  particles  have  been  thoroughly  investigated  both 
experimentally  and  theoretically  [1,  2,  3,  4,  5]  as  they  can  give  useful  information  e. 
g.  on  the  nature  and  on  the  size  of  this  kind  of  scatterers.  The  resonances  may  give 
useful  information  even  for  non  spherical  particles  because  the  resonant  behavior  is 
expected  to  be  strongly  dependent  on  morphology  as  well  as  on  the  orientation  of  the 
particles  with  respect  to  the  incident  field. [6,  7]  For  instance,  the  literature  reports 
studies  intended  to  investigate  the  suppression  of  some  of  the  resonances  in  the  spectrum 
of  axially  symmetric  particles  for  particular  choices  of  the  direction  of  incidence  and  of 
the  polarization  of  the  incident  wave. [8,  9] 

Among  the  methods  that  are  applicable  to  calculate  the  extinction  spectrum  of 
particles,  the  transition  matrix  approach  [10]  plays  a  prominent  role.  Indeed,  the 
related  formalism,  together  with  the  theory  of  the  vector  multipole  fields, [11]  allows 
for  the  separation  of  the  effects  that  are  due  to  the  morphology  of  the  particles  from  the 
effects  that  are  due  to  their  orientation  with  respect  to  the  incident  field.  As  a  result  it 
becomes  an  easy  matter  to  perform  analytical  angular  averages  so  as  to  get  information 
on  the  optical  behavior  of  low-density  disj^ersions  with  known,  or  assumed,  distribution 
of  the  orientations  of  the  particles.  [12]  Such  averaging  technique  has  been  applied 
to  calculate  the  extinction  spectrum  from  dispersions  both  of  aggregates  of  spherical 
scatterers  and  of  homogeneous  spheres  containing  spherical  inclusions. [12,  13,  14]  As 
individual  scatterers  the  latter  are  particularly  interesting  because  they  may  well  appear 
as  spherically  symmetric  but  their  calculated  spectrum  displays  a  noticeable  dependence 
on  the  polarization  and  on  the  direction  of  the  incident  field.  Even  in  the  simple  case 
of  spheres  containing  a  single  inclusion,  a  calculation  of  the  optical  resonances  may  be 
useful  to  put  into  evidence  the  kind  of  the  changes  undergone  by  the  spectrum  when  the 
inclusion  occupies  more  and  more  eccentric  positions.  More  precisely,  such  calculations 
may  help  to  assess  whether  an  experimental  resonance  spectrum  may  be  effective  to 
probe  the  internal  asymmetry  that  is  due  to  the  eccentricity  of  the  inclusion  as  well  as 
to  give  indication  on  the  fraction  of  the  scattering  particles  that  are  likely  oriented. 

With  this  in  mind,  the  present  paper  is  devoted  to  the  calculation  of  the  extinction 
spectrum  of  a  homogeneous  sphere  of  MgO  containing  a  spherical  inclusion  of  Mg  for 
a  few  different  positions  from  the  centered  to  the  most  eccentric  one.  Our  choice  of 
the  material  is  not  accidental  as  metal  particles  eccentrically  coated  by  their  oxyde  are 
interesting  e.  g.  in  astrophysics  as  their  presence  in  the  interstellar  dust  is  invoked  to 
explain  the  observed  dependence  on  the  polarization  of  the  galactic  extinction. [15]  The 
electromagnetic  scattering  properties  of  this  kind  of  particles  are  calculated  by  resorting 
to  the  technique  of  Refs.  [12,  13,  14].  Both  the  incident  and  the  scattered  field  as  well  as 
the  field  within  the  particle  are  described  through  the  appropriate  expansions  in  terms 
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of  vector  multipole  fields  whose  amplitudes  are  determined  by  the  boundary  conditions 
across  the  surface  of  each  component  of  the  scatterer.  In  fact,  this  procedure  enables 
us  to  calculate  the  transition  matrix  as  a  function  of  the  position  of  the  inclusion,  then 
the  scattered  field  for  given  polarization  and  direction  of  incidence  of  the  incoming  field. 
Furthermore,  by  resorting  to  the  well  defined  transformation  properties  of  the  multipole 
fields  under  rotation,  we  perform  an  average  over  the  angular  distribution  of  the  particles 
to  get  information  on  the  resonant  behavior  of  a  dispersion  of  particles. 

In  Section  2  we  introduce  the  multipole  expansion  of  the  electromagnetic  field  and 
define  the  elements  of  the  scattering  amplitude  in  terms  of  the  transition  matrix  elements 
and  of  the  multipole  amplitudes  of  the  incident  field. 

In  Section  3  we  present  the  results  of  our  calculation  of  the  extinction  efficiency  of 
a  dielectric  sphere  of  MgO  containing  a  spherical  inclusion  of  Mg. 

A  few  concluding  remarks  will  be  drawn  in  Section  4. 


2.  Transition  matrix  and  scattering  amplitude 

We  partition  the  space  into  the  three  regions  sketched  in  Fig.  1:  the  external  region,  that 
is  assumed  to  be  filled  by  a  homogeneous,  non-dispersive,  non-absorbing  medium  with 
refractive  index  n  (typically  the  vacuum);  the  interstitial  region,  centered  at  the  origin 
and  of  radius  that  is  filled  with  a  homogeneous,  possibly  dispersive  and  absorbing 
medium  with  refractive  index  no]  the  region  within  the  inclusion,  centered  at  Ri  and  of 
radius  /?i,  characterized  by  a  refractive  index  that  may  be  dispersive  and  absorbing. 
For  the  sake  of  simplicity,  the  theory  in  this  section  will  refer  to  homogeneous  inclusions 
only,  on  account  that  the  resulting  formulas  are  easily  extended  to  the  case  of  radial 
nonhomogeneity.  [17] 

We  assume  that  all  the  fields  depend  on  time  through  the  factor  and  define 

the  propagation  constants 

K  =  kn^  Kq  =  kno^  Ki  =  kni , 

in  each  of  the  regions  mentioned  above,  respectively,  with  k  =  u;/c.  The  incident  field 
is  assumed  to  be  the  plane  wave  of  amplitude  Eq 

Elr)  —  EoXXir)  , 

where  K/  =  A^k/  is  the  incident  wavevector;  is  a  unit  polarization  vector  whose 
index  rj  indicates  whether  the  field  is  parallel  (?/  =  1)  or  perpendicular  {rj  =  2)  to 
the  plane  of  scattering,  i.  e.  to  the  plane  that  contains  both  K/  and  the  direction  of 
observation.  The  unit  vectors  Ujrj  are  so  oriented  that 

U/i  X  U/2  ^  k/. 
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Since  our  approach  makes  extensive  use  of  the  multipole  expansion  of  the 
electromagnetic  field,  let  us  define  the  multipole  fields  as 

where  the  functions  X^Af  vector  spherical  harmonics  [18]  and,  on  account  of  the 
chosen  time  dependence  of  the  fields,  the  quantities  are  spherical  Hankel  functions 
of  the  first  kind;  the  superscript  p  is  a  parity  index  that  distinguishes  the  magnetic 
multipoles  (p  =  1)  from  the  electric  ones  (p  =  2).  We  also  need  to  define  the  multipole 
fields  that  are  identical  to  the  fields  Eq.  (1),  except  for  the  substitution  of 

the  spherical  Bessel  functions,  in  place  of  the  Hankel  functions,  hi. 

The  multipole  expansion  of  the  interstitial  field  as  well  as  of  the  field  within  the 
inclusion  are  fully  described  in  ref.  [13].  Here  we  report  only  the  expansion  of  the  field 
in  the  external  region  where  it  is  the  superposition  of  the  incident  and  of  the  scattered 
field  and  is,  accordingl}^,  expanded  as 

E„  =  E5,  +  E„  =  E  {/iaH«{r.  K)  +  (A/..  k;)jS(r.  /O}  ■ 

plm 

The  multipole  amplitudes  of  the  field  scattered  by  the  whole  object,  are  expected 

to  depend  both  on  the  polarization  and  on  the  direction  of  incidence  a.s  well  as  on 
the  position  of  the  inclusion  within  the  external  sphere.  The  quantities  Wi^  are  the 
multipole  amplitudes  of  the  incident  field  that  for  a  plane  wave  with  polarization  along 
the  unit  vector  u  are  defined  as  [9,  12] 


h1m(u,  k)  =  47r2^u  •  Xz,A/(k),  ^  ^2) 

H^Ak^^k)  =  47ri^+^(k  x  u)  •  XLA^(k). 

The  amplitudes  are  the  main  unknowns  of  the  problem  and  are  determined  by 

boundary  conditions  at  the  surface  of  the  inclusion  as  well  as  at  the  surface  of  the 
external  sphere.  The  procedure  for  calculating  the  A-amplitudes  is  fully  described  in 
Ref.  [13].  For  our  present  purposes  we  only  need  to  recall  that  the  W  and  the  A- 
amplitudes  are  related  by  the  equation 


(3) 


where  the  quantities  are  the  elements  of  the  so  called  transition  matidx  S, 

that  accounts  for  the  morphology  (geometry  and  scattering  power)  of  the  particle.  [10] 
Equation  (3)  shows  that  the  resonances  occur  when,  with  varying  frequency  of  the 
incident  wave,  at  least  one  of  the  elements  of  S  becomes  much  larger  than  any  other 
element. [19]  As  a  result  some  of  the  amplitudes  A^^^^  become  large  thus  yielding  a 
resonance  peak  that,  for  a  given  polarization,  ?/,  can  be  classified  by  the  labels  p,  /,  and 
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m.  Nevertheless,  however  large  the  elements  of  S  may  be,  no  resonance  may  occur  if 
the  appropriate  amplitudes  in  the  expression  of  the  incident  field  are  vanishing. 

The  explicit  expression  of  the  transition  matrix  for  a  sphere  containing  an  eccentric 
inclusion  can  be  easily  obtained  by  solving  formally  Eqs.  (6)  and  (7)  of  Ref.  [13].  The 
result  is 


S  —  (Mo  —  Mvi/lo^iRili^o)(Ro^ 

The  elements  of  the  matrices  in  the  right  hand  side  of  Eq.  (5)  are  explicitly  given  in 
Ref.  [13].  All  these  matrices  are  diagonal  except  for  Iq^i  and  h^o  that  account  for  the 
eccentricit}^  of  the  inclusion:  As  a  result  S  is  a  non-diagonal  matrix.  Nevertheless,  in 
order  to  get  a.  simpler  exj^ression  lor  the  transition  matrix  we  assume  that  the  diameter 
along  which  the  center  of  the  inclusion  lies  coincides  with  the  z  axis.  With  this  choice 
of  the  geometry,  in  fact,  the  cylindrical  symmetiy  around  the  2:  axis  makes  the  elements 
of  S  to  vanish  for  m  ^  m\  Accordingl}^,  they  will  be  hereafter  indicated  as 
Eq.  (3)  takes  on  the  simpler  form 


p'l* 


It  may  be  useful  to  remark  that  when  Ri  — ^  0,  i.  e.  when  the  inclusion  is  centered, 
the  scatterer  becomes  spherically  S3m'imetric  and  its  transition  matrix  becomes  diagonal 
with  elements  independent  of  in.  In  fact,  it  may  be  easily  verified  that  the  limiting 
expression  of  the  matrix  S  coincides  with  the  expression  that  is  appropriate  for  a  sphere 
containing  a  concentric  inclusion.  [20] 

As  usual,  we  describe  the  scattered  field  in  the  far  zone  through  the  normalized 
scattering  amplitude  that  is  defined  by  the  equation  [18] 

The  arguments  of  f,,  recall  that,  in  general,  the  scattering  amplitude  depends  both  on 
the  incident  and  on  the  observed  wavevector  as  well  as  on  the  polarization,  ?/,  of  the 
incident  field.  In  terms  of  the  amplitudes  of  the  scattered  field  the  expression  of  is  [13] 

f.  =  i  E(-0'”l>lSi.X,„(k5)  +  iA'-Z^s  X  X,„(t5)|,  (6) 

Im 

where  £5  denotes  the  direction  of  observation.  Of  course,  the  scattered  field  is  neither 
parallel  nor  perpendicular  to  the  plane  of  scattering:  It  can  be  decomposed,  however, 
into  its  components  that  are  parallel  and  perpendicular  to  this  plane  by  projecting 
on  the  pair  of  unit  vectors  Usp  such  that  Us2  =  fi/2  ^nd 


U51  X  U52  ~  ks- 
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These  vectors  form  the  polarization  basis  that  is  appropriate  for  the  scattered  wave. 
The  components  of  along  these  vectors  turn  out  to  be  given  by  the  equation 


fr)r}'  —  ^r)  *  ^Sr)' 

'  I E  E  E  (o,„  k,) 


iirnk 


pi  p*V  rn 


“  ^vT]'h 


(7) 


that  is  easily  obtained  from  Eq.  (6)  with  the  help  of  Eqs.  (2)  and  (5).  The  2x2  matrix  of 
elements  7,^/  has  been  used  by  van  de  Hulst  [21]  to  illustrate  the  symmetry  properties  of 
a  general  scattering  process.  In  particular,  these  symmetry  properties,  when  considered 
together  with  the  cylindrical  symmetry  of  the  scatterer,  imply  that  —  0  when  ?/  ^ 
viz.  the  scafterers  we  are  dealing  with  here  never  give  rise  to  cross-polaxization  effects. 
Equation  (7)  shows  that  depends  not  only  on  the  properties  of  the  transition 
matrix  S  but  also  on  the  proi:)erties  of  the  M'^-amplitudes,  In  order  to  simplify  our 
considerations  we  assume  that  the  plane  of  scattering  coincide  with  the  x-z  plane;  this 
choice,  however  does  not  imply  any  loss  of  generality  because  choosing  any  other  plane 
through  the  2:  axis  will  onl}^  introduce  an  inessential  phase  factor.  Now,  with  the  help 
of  the  definition  of  the  M^-amplitudes,  Eq.  (2),  it  is  an  easy  matter  to  prove  that 

k;)  =  k/),  (8) 


whereas  the  transition  matrix  elements  have  the  property 


^(p^)  _  ^_^p-p' 


c(PiP0 


(9) 


that  is  inferred  from  their  explicit  expression.  Equations  (8)  and  (9)  ensure  that  the 
scattered  field  have  the  correct  transformation  properties:  In  fact,  the  true- vector  nature 
of  the  electric  field  requires  that  it  remain  unchanged  under  reflection  in  the  plane  of 
scattering  for  parallel  polarization  (7/  =  1),  whereas  it  must  change  its  sign  for  orthogonal 
polarization  (7/  =  2).  An  useful  implication  of  Eqs.  (8)  and  (9)  is  that  the  resonances 
for  a  given  77,  p,  and  /  that  are  associated  to  ±77z  are  degenerate,  i.  e.  they  must  occur 
at  the  same  frequency. 

We  finally  emphasize  that  Eq.  (7)  presents  all  the  advantages  of  an  approach  based 
on  the  transition  matrix  and  on  the  multipole  expansion  of  the  fields.  Equation  (7), 
in  fact,  describes  the  scattering  from  a  particle  of  given  orientation.  Nevertheless,  by 
resorting  to  the  transformation  properties  of  the  PE-amplitudes  under  rotation,  it  is 
an  easy  matter  to  perform  analiticaJ  averages  over  the  orientation, [13]  The  averaging 
procedure,  that  has  been  fully  discussed  elsewhere, [12]  yields  information  on  the 
extinction  of  a  low  density  dispersion  of  particles,  and  thus  produces  results  that  can 
be  directly  compared  with  the  experimental  data. 
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3.  Results  and  discussion 

The  theory  in  Section  2  has  been  applied  to  the  calculation  of  the  extinction  spectrum 
of  spheres  of  MgO  with  a  radius  po  —  25  nm  containing  a  spherical  inclusion  of  Mg  with 
a  radius  of  pi  =  16.6  nm.  In  other  words  we  are  considering  a  metal  sphere  coated  by  its 
oxide,  po  ^nd  pi  were  chosen  so  as  to  lie  within  the  estimated  ranges  for  the  particles  in 
the  interstellar  medium. [15]  Furthermore,  our  choice  implies  that  the  particles,  although 
small,  can  have  resonances  also  for  /  >  1.  The  dielectric  function  of  Mg  was  assumed  to 
be  of  the  free-electron  Drude  form  [22] 


uj{uj  +  iVu.,y 

where  =  1.6  •  lO^^sec"^  and  T  —  0.01.  In  turn,  the  dielectric  function  of  MgO  was 
assumed  of  the  damped  oscillator  form  [22] 

€  -  Coo  __  _ _ 

Co  “Coo  ’ 

with  UJT  =  7.5  •  10^^sec-\  F  =  0.27,  6o  9.8  and  Coo  =  2.95.  The  results  of  our 
calculations  are  presented  by  reporting  the  extinction  efficiency 
4  >.  .  . 

Po 

as  a  function  of  the  ratio  i/  =  u/up.  In  the  preceding  equation  £5  =  k/)  is  the 

forward-scattering  amplitude  for  incidence  along  k/  =  (??/,(p/).  Of  course  when  the 
scatterers  are  sphericall}^  S3anmetric  Qy^  is  actually  independent  of  the  polarization  and 
the  index  t]  will,  accordingly,  be  dropped. 

In  Fig.  2  we  report  Q  for  a  sphere  of  solid  Mg  with  radius  pi.  In  order  to  check 
the  convergency  of  our  calculations  and  to  associate  the  computed  resonances  to  the 
appropriate  value  of  /,  the  extinction  efficiency  is  reported  for  /m  =  1  Im  —  2,  Im 
being  the  maximum  value  of  /  included  in  the  sums  in  Eq.  (7).  We  notice  that  taking 
>  2  does  not  cause  visible  changes  on  the  scale  of  the  figure  nor  yields  the  occurrence 
of  any  higher  order  resonance.  However,  in  order  to  achieve  a  numerical  accuracy  to 
four  significant  digits  all  our  calculations  were  actually  performed  with  Im  =  5.  The 
sphere  of  Mg  presents  two  resonances:  the  lowest-frequency  peak  is  associated  with 
/  =  1  whereas  the  second  peak  is  associated  with  I  =  2.  An  analysis  of  the  elements  of 
the  transition  matrix  suggests  that  both  peaks  are  due  to  electric  resonances  and,  due 
to  the  independence  of  in  of  the  elements  of  the  transition  matrix  for  a  sphere, [21]  the 
value  of  m  need  not  be  specified.  We  do  not  report  Q  for  the  sphere  of  MgO  of  radius 
po  because  this  scatterer  does  not  present  any  resonance  in  the  range  of  interest. 

Our  calculations  for  the  actual  scatterer  were  performed  with  the  geometry  described 
in  Section  2,  for  several  positions  of  the  inclusion  from  the  centered  position  to  the  most 
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eccentric  one,  so  as  to  follow  the  evolution  of  the  resonances  a.s  the  inclusion  becomes 
more  and  more  eccentric.  As  one  should  expect,  we  found  that  the  effects  that  are  due  to 
the  eccentricity  achieve  the  maximum  of  evidence  when  the  inclusion  becomes  tangent 
to  the  surface  of  the  external  sphere  so  that  in  the  following  figures  only  the  results 
for  this  extreme  case  are  reported.  Accordingly,  Fig.  3  shows,  for  incidence  along  the 
z  axis,  Qrj  for  the  actual  scatterer  with  the  inclusion  at  maximum  eccentricity  together 
with  the  partial  efficiencies 
4 
Po 

with  /  =  1,2.  Of  course,  due  to  the  cylindrical  symmetry  of  the  scatterer,  both  the 
total  and  the  partial  extinction  efficiencies  are  independent  of  ?/.  Figure  3  also  reports 
Q  for  the  centered  inclusion.  In  the  latter  ca,se  a  single  electric  resonance  occurs  that 
is  associated  to  /  —  1  whereas  the  resonance  that  for  a  sphere  of  solid  Mg  is  associated 
to  I  —  2  turns  out  to  be  so  damped  as  to  reduce  to  a  small  shoulder.  This  behavior 
is  easily  understood  in  the  framework  of  the  theory  of  Johnson:  [23]  On  account  of 
the  analogy  between  the  radial  Schrodinger  equation  for  spherically  symmetric  s^^stems 
and  the  radial  Helmholz  ecjuation  for  spherical  scatterers,  this  author  interpretates  the 
electromagnetic  resonances  as  quasi-bound  states  of  the  radiation  within  the  effective 
potential  well 

The  presence  of  a  coating  may  modify  the  radial  shape  of  the  refractive  index  72(r)  to 
such  an  extent  that  the  potential  well  for  the  appropriate  value  of  /  (in  the  present  case 
for  1  =  2)  may  not  occur  at  all,  thus  preventing  the  occurrence  of  the  corresponding 
resonance,  or  may  become  so  shallow  as  to  produce  the  damping  that  is  observed  in 
Fig.  3.  Considerations  of  this  kind  do  not  apply  when  the  inclusion  is  off  center  because 
of  the  lack  of  spherical  S3anmetry.  Indeed,  Fig.  3  shows  that  when  the  inclusion  is  off 
center  the  spectrum  complicates  as  four  resonances  now  occur.  The  one  at  the  lowest 
frequency  is  associated  to  /  =  I  and  is  the  evolution  of  the  main  peak  that  occurs 
when  the  inclusion  is  centered.  The  three  further  resonances  can  be  attributed  to  the 
multiple  scattering  processes  between  the  inclusion  and  the  external  sphere.  [13]  Due 
to  the  off-center  position  of  the  inclusion  these  processes  are  asymmetric  and  thus  yield 
the  observed  multiplication  of  peaks.  The  behavior  of  the  quantities  qi  clearly  shows 
that  all  these  peaks  arise  from  the  superposition  of  contributions  for  /  =  I  and  /  =  2. 
The  occurence  of  this  superposition  calls  for  a  word  of  caution  about  the  classification 
of  the  resonances  into  electric  and  magnetic  as  well  as  about  the  order  of  the  resonating 
multipoles.  In  a  previous  paper  [24]  we  noticed  that  this  classification  depends  on  the 
origin  to  which  the  multipole  fields  are  referred.  For  a  sphere  with  a  centered  inclusion 
there  exist  a  natural  origin,  namel}^  the  center  of  the  external  sphere  that  coincides 
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with  the  center  of  symmetry  of  the  whole  scatterer;  of  course,  when  the  inclusion  is 
eccentric  such  a  center  of  symmetry  does  not  exist.  Now,  a  multipole  field  that  is 
purely  electric  or  magnetic  when  referred  to  a  given  origin  becomes  a  superposition  of 
both  electric  and  magnetic  multipole  helds,  in  general  of  all  the  multiplicities,  when 
referred  to  a  different  origin.  [25]  On  the  other  hand,  in  the  framework  of  our  present 
approach,  the  addition  theorem  of  Ref.  [25],  gives  rise  to  the  non-diagonal  matrices  Iq^i 
and  li^o  ill  the  expression  of  the  transition  matidx.  As  these  matrices  account  for  the 
multiple  scattering  processes  within  the  components  of  the  scatterer,  we  feel  authorized 
to  state  that  the  high  frequency  peaks  are  due  to  multiple  scattering  processes  and  that 
their  electric  or  magnetic  nature  cannot  be  stated  in  an  absolute  sense  but  only  with 
reference  to  a  chosen  origin  that,  in  the  present  case,  coincides  with  the  center  of  the 
external  sphere.  With  tliis  exception  in  mind,  we  can  safely  state,  however,  that  all 
the  resonances  in  Fig.  3  belong  to  |?7i|  =  1  becanse,  for  any  ?/,  p,  and  /,  the  multipole 
amplitudes  of  a  plane  wave  that  propagates  along  the  2:  axis  do  not  vanish  for 

m  =  ±1  onl}^  Furthermore,  according  to  our  remark  in  Section  2,  Eqs.  (8)  and  (9)  imply 
that  the  resonances  associated  to  m  and  to  —m  must  occur  at  the  same  frequency.  We 
finally  notice  that  the  same  Siiectrum  is  obtained  by  changing  the  direction  of  incidence 
from  di  =  0°  to  iSj  =  180°:  This  is  an  expected  result  that  stems  from  the  symmetry 
properties  of  fr^rj  for  forward  scattering.  [21] 

The  results  that  we  presented  so  far  show  the  expected  independence  of  the 
polarization  either  because  the  scatterer  is  spherically  symmetric.  Fig.  2,  or  because 
the  incidence  is  along  the  C3dindrical  axis.  Fig.  3.  In  this  respect  let  us  recall  that, 
according  to  Eq.  (7),  the  dependence  of  the  scattering  amplitude  on  the  polarization  as 
well  as  on  k/  and  on  ks  is  entirel}^  due  to  the  amplitudes  that  are  independent  of 
frequency,  however.  On  the  contrary,  the  possibility  of  producing  resonances  is  contained 
in  the  frequenc3^-dependence  of  the  elements  of  the  transition  matrix,  that,  in  turn  do 
not  depend  either  on  77  or  on  k/  or  on  k5.  With  this  in  mind,  let  us  go  to  discuss  the 
results  of  Fig.  4  where  the  dependence  of  the  spectrum  on  the  choice  of  the  polarization 
shows  to  its  full  extent.  In  Fig.  4,  indeed,  we  superpose  to  the  spectrum  for  a  sphere 
with  a  centered  inclusion  (the  same  that  we  reported  in  Fig.  3)  the  extinction  efficiency 
for  the  case  of  the  inclusion  at  maximum  eccentricity  for  polarization  both  parallel 
and  perpendicular  to  the  plane  of  scattering;  the  direction  of  incidence  is  individuated 
by  the  polar  angles  7?/  =  52°,  (pj  =  0°.  This  choice  of  the  incidence  does  not  bear 
any  particular  significance  thus  ensuring  the  true  generality  of  our  results.  In  the 
same  figure  we  also  report  the  result  of  the  angular  average  over  the  particles  with  the 
eccentric  inclusion  on  the  assumption  that  they  are  randomly  distributed  in  orientation; 
of  course,  the  latter  specti-um,  that  gives  information  on  the  extinction  from  a  dispersion 
of  identical  scatterers,  is  independent  of  the  polarization.  We  first  notice  that,  except 
for  the  leftmost  peak,  all  I'esonances  for  the  case  of  an  eccentric  inclusion  arise  from  a 
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superposition  of  contributions  froin  /  =  1  and  for  I  =  2;  this  is  shown  the  values 
of  the  partial  efficiencies  cind  tfia.t  are  not  reported,  however.  Therefore  our 
previous  considerations  both  on  tlie  origin  and  on  the  classification  of  the  resonances 
apply.  Comparison  of  the  spectra  lor  the  individual  particles  shows  that  the  resonance  at 
i/j  =  0.535  occurs  only  in  j)erpendicular  polarization  whereas  the  resonance  at  //2  =  0.545 
occurs  only  in  parallel  polarization.  The  separation  U2  ~  although  small  on  the  scale 
of  the  figure,  is  quite  noticeable  in  wavelength.  The  peak  at  i^i  is  the  same  that  appears 
in  Fig.  3  for  i?/  =  0°.  In  fact  this  resonance  mainly  stems  from  the  elements 
that,  for  =  0°,  are  excited  for  both  polarizations  as  the  relevant  amplitudes  kFi,2i 
are  independent  of  the  choice  of  ?/.  On  the  contrary,  for  general  direction  ol  incidence 
1Fi^2i  unchanged  for  i]  —  2  whereas  they  are  smaller  lor  7/  =  1.  The  peak  at  1^2 
originates  instead  from  the  leading  element  excited  for  p^irallel  polarization 

as  the  exciting  amplitude  provided  that  the  incidence  is  not  along  the  2:  axis,  does 

not  vanish  for  ?/  =  1  only. 

When  the  particles  are  randomly  oriented,  the  averaged  spectrum  shows  that  both 
peaks  at  //i  a.nd  U2  <^re  present.  This  is  an  expected  result  that  we  already  noticed  just  in 
connection  with  the  optical  properties  of  sj^heres  containing  an  eccentric  inclusion. [13] 
Averaging  over  the  orientations,  indeed,  implies  summing  the  contributions  to  the 
extinction  efficiency  from  scatterers  whose  symmetry  axis  does  not  lie,  in  general,  in  the 
plane  of  scattering.  Now,  the  resj^onse  of  a  scatterer  with  its  s3umTietry  axis  not  lying  in 
the  scattering  plane  is  easil}^  recognized  as  identical  to  the  response  of  a  scatterer  with 
its  axis  in  that  plane  when  this  latter  object  is  excited  by  a  wave  with  an  appropriate 
state  of  polarization  that,  in  general,  is  neither  parallel  nor  orthogonal:  Therefore  the 
above  mentioned  response  must  be  a  linear  combination  of  the  responses  for  parallel 
and  orthogonal  polarization.  Since  the  extinction  efficiency  of  the  individual  scatterers 
displays  a  peak  whose  frequency  changes  with  the  polarization  we  get  the  spectrum  for 
random  orientation  that  is  reported  in  Fig,  4.  Ultimately,  at  least  one  of  the  resonances 
of  the  individual  scatterers  undergoes  evident  qualitative  changes  as  a  function  of  the 
polarization.  Thus  the  appearence  of  a  dependence  on  the  polarization  in  the  resonance 
spectrum  from  a  dispersion  of  such  particles  may  provide  an  indication  that  a  noticeable 
fraction  of  them  is  likel}^  oriented. 

4,  Conclusion 

The  results  that  we  presented  in  section  3  suggest  that  a  careful  anal^^sis  of  the  extinction 
spectrum  from  spheres  containing  a  spherical  inclusion  may  help  establishing  whether 
the  inclusion  is  concentric  or  eccentric.  In  the  latter  case,  in  fact,  the  eccentricit}^  yields 
a  multiplicit}^  of  resonances  that  are  quite  absent  in  the  spectrum  from  spheres  with  a 
centered  inclusion.  In  the  case  that  the  inclusion  is  eccentric,  measurements  performed 
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with  polarized  light  may  give  indication  on  the  fraction  of  particles  that  are  oriented 
alike.  Getting  this  kind  of  information  can  be  of  importance  e.  g.  in  astrophysics.  Indeed, 
the  galactic  extinction  shows  an  evident  dependence  on  the  polarization  for  certain  lines 
of  sight.  This  fact  has  been  interpreted  as  due  to  the  presence  of  anisotropic  particles 
of  the  kind  we  dealt  with  in  this  paper  which  are  partially  oriented  by  some  magnetic 
field.  Except  for  the  cause  of  the  allignment,  that  does  not  concern  us  here,  the  results 
in  Section  3  support  the  reliability  of  this  interpretation. 
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Figure  captions 


Figure  1.  Sketch  of  the  geometry  for  a  sphere  with  an  eccentric  inclusion.  In  actual 
calculations  the  center  of  the  inclusion  was  assumed  to  lie  on  the  z  axis. 


Figure  2.  Extinction  efficiency  for  a  sphere  of  solid  Mg  with  radius  p\  —  16.6  nm. 
The  continuous  curve  was  obtained  with  Im  ~  I  and  the  dashed  curve  with  =  2. 


Figure  3.  Extinction  efficiency  for  the  actual  scatterer  with  the  inclusion  centered 
(continuous  curve  with  dots)  and  with  the  inclusion  at  maximum  eccentricity 
(continuous  curv'^e).  The  dashed  and  the  dotted  curves  report  the  partial  efficiencies 
for  /  —  1  and  /  —  2,  respectively,  for  the  case  of  maximum  eccentricity.  The  incidence 
is  along  the  z  axis  {'df  =  0°,  =  0°). 


Figure  4.  Extinction  efficiency  for  the  actual  scatterer  with  the  inclusion  at  maximum 
eccentricity  for  polarization  parallel  (dashed  curve)  and  perpendicular  (dotted  curve) 
to  the  plane  of  scattering.  The  incidence  is  individuated  by  the  polar  angles  i?/  — 
52°,  <pj  =  0°.  The  continuous  curve  reports  the  extinction  efficiency  of  a  dispersion  of 
scatterers  with  random  distribution  of  the  orientations.  A  detail  of  the  rightmost  part 
of  the  spectrum  is  shown  in  the  inset.  For  comparison  the  efficiency  for  the  case  of  the 
centered  inclusion  (continuous  curve  with  dots)  is  also  reported. 
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Figure  4.  Extinction  efficiency  for  the  actual  scatterer  with  the  inclusion  at  maximum 
eccentricity  for  polarization  parallel  (dashed  curve)  and  perpendicular  (dotted  curve) 
to  the  plane  of  scattering.  The  incidence  is  individuated  by  the  polar  angles  dj  = 
52®,  (fj  =0®.  The  continuous  curve  reports  the  extinction  efficiency  of  a  dispersion  of 
scatterers  with  random  distribution  of  the  orientations.  A  detail  of  the  rightmost  part 
of  the  spectrum  is  shown  in  the  inset.  For  comparison  the  efficiency  for  the  case  of  the 
centered  inclusion  (continuous  curve  with  dots)  is  also  reported. 


